Two-dimensional electron scattering in regions of nonuniform spin-orbit coupling 
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We present a theoretical study of elastic spin-dependent electron scattering caused by a nonuni- 
form Rashba spin-orbit coupling strength. Using the spin-generalized method of partial waves the 
scattering amplitude is exactly derived for the case of a circular shape of scattering region. We found 
that the polarization of the scattered waves are strongly anisotropic functions of the scattering an- 
gle. This feature can be utilized to design a good all-electric spin-polarizer. General properties of 
the scattering process are also investigated in the high and low energy limits. 
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The presence of the spin-orbit interaction (SOI) de- 
stroys the spin-rotational symmetry, therefore the prop- 
erties of the scattering get influenced by the spin state of 
the incident particle^. The SOI may lead to the asym- 
metry of the differential scattering cross section (skew- 
scattering), and it may affect the polarization vector of 
the incident hea.w^'^. In nuclear physics, this prop- 
erty has been utilized to generate spin-polarized neutrons 
from an unpolarized beam by scattering it on a zero- 
spin nucleus^. In low-dimensional semiconductors, sig- 
nificant spin-splitting in the absence of a magnetic field 
is observed, which is mostly attributed to the SOI of 
the Rashba type arising from the structural inversion 
asymmetry of the hosting heterostructure^. Promising 
spin transistor application has been proposed exploiting 
the tunability of the strength of the Rashba coupling by 
an external electrostatic field^, and initiated an inten- 
sive research in the field of spintronics^. To generate 
spin-polarized electron beams, which is a fundamental 
problem in spintronics, several mechanisms have been 
proposed". 

Recently, all-electrical (without externally applied 
magnetic fields) spin-polarizer devices have been 
suggested^ utilizing the properly designed spatial modu- 
lation of the Rashba SOI strength a. The spatial vari- 
ation of the Rashba strength, which is proportional to 
the magnitude of the electric field applied perpendicular 
to the two-dimensional electron gas (2DES) systems, can 
be achieved by small biased electrodes on the top of the 
heterostructurc. 

In this work we show that the polarization of the elas- 
tically scattered wave caused by a nonuniform Rashba 
strength becomes strongly anisotropic, ie depends on the 
angle (called scattering angle) between the directions of 
the incident electron beam and that of the scattering 
wave. We also demonstrate that (i) the differential scat- 
tering cross section has a skew-scattering featureiiSi^ for 
polarized incoming electron beams, (ii) using experimen- 
tally relevant material parametersifl an almost full po- 
larization of unpolarized incident electron beams can be 
observed in a narrow window of scattering angles. More- 
over, our analytical calculations allows us to derive uni- 
versal properties for the scattering amplitude and the 



polarization P*"^ in the high and low energy limits. 

To this end, we consider a system in which the Rashba 
SOI strength a{r) varies on the plane of the 2DES as 
a{r) = aiQ{a — |r|) -f Q;20(|r| — a), where a is the radius 
of the scattering center, O is the Heaviside function and 
r = {x,y) defines the position on the plane (see Fig. 
The Hamiltonian of the system in the one-band effective- 




FIG. 1: (Color online) The plane wave of an incident electron 
with wave vector k gets scattered by a scattering region de- 
fined by a nonuniform Rashba coupling strength a(r) (see the 
text), and its portion propagates along the direction given by 
the scattering angle ip. 
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where p — {px,Py) is the momentum operator, m* is the 
effective mass of the electron, , cr^ are the Pauli matri- 
ces. Note that this Hamiltonian is symmetrized to make 
it Hermitian. For simplicity, we set a2 = 0, and ai is a 
constant value throughout this paper (our analysis can 
easily be extended to the case a2 7^ 0). Recently, similar 
scattering problems have been studied with uniform a(r) 
and an electrostatic potential varying in the plane of the 
2DES"i2. 

The spin-density matrijsSiSi is proved to be useful for 
treating the coupled spin-charge quantum transport in 
spintronic devices Employing this formalism allows 
us to derive explicit formulas for the differential and the 
total scattering cross section, and the polarization P^"^ of 
the scattered wave in terms of the polarization P'"^ of 
the incident electron beam. 

As can be seen below, the physics of the scattering of 
electrons shown in Fig. ^ is fundamentally different from 
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the Mott scattering originating from the spin-dependent 
scattering potentialiiSiM. and the spin-dependent two- 
dimensional electron scattering from quantum dots and 
antidots studied in Ref. Q. In these scattering prob- 
lems the polarization P*"^ of the scattered wave is al- 
ways perpendicular to the scattering plane for unpolar- 
ized (P™'^ = 0) incident particles of spin 1/2, while in 
our scattering problem this is not the case in general. 

First, we briefly outline the spin-density matrix ap- 
proach of our spin-dependent scattering problem. The 
incident plane wave with wave number k has the form 
tp"^^[r) = e^^'^ I7) , where I7) denotes the spin state. In 
two dimensions the form of the scattered wave asymptot- 
ically far from the scattering center is 



(2) 



where w = 2Re(uQu) + i{u x u*) and the components of 
the matrix Al: 



M^J = {\uof ~\uf)S,j+2Re{u*Uj) 
3 

+ 2y^gij7clm (uQUk), 



(6b) 



with i,j = 1, 2, 3, and Sij and Sijk denote the Kronecker 
delta and the Levi-Civita symbol, respectively. Here 
Im(-) stands for the imaginary part of the argument. The 
components of Al are real numbers. 

The spin-dependent form of the optical theorem can 
be derived by considering the scattering time-evolution 
of Gaussian wave packets, and we obtain: 



CTtot 



[«o(0) + u(0).P'-]}, (7) 



where t{ip) is the scattering amplitude (2x2 matrix in 
spin space), and depends on the scattering angle ip and k. 
The scattering amplitude t{ip) can be expanded in terms 
of the unit matrix (Tq and the vector cr = (cri, tT2, fs) 
formed from the three Pauli matrices (for convenience, 



we use (Ti = a-x, cr2 
Pauli matrices): 



cTj, , (T3 = az notations for the three 



(3) 



fe=0 



where ^-nd u = (mi, 1*2,^3) can only be obtained by 
solving the Schrodinger equation for the scattering prob- 
lem. 

The spin-density matrix p™'^ = i (ctq + c ■ P™'^) of the 
incident wave for a given polarization P'"'^ = (cr)inc — 
Tr [p™'^ a) is related to the spin-density matrix p^'^ of the 
scattered wave as^ 



f^incft 

Tr(fpi"<=ft)' 



(4) 



where Tr denotes the trace in the spin states. Then using 
(O the differential scattering cross section reads as 



da 

dip 



= Tr (fp'^'^ft) = c + V • P'"^ where 



(5a) 



k=0 



\ukf and V = 2Re(MQu) - z(u X u*). (5b) 



Here Re(-) denotes the real part of the argument, and 
the star stands for the complex conjugation. Note that, 
in general, u x u* is not zero but it is always a purely 
imaginary vector. 

Similarly, we found for the polarization vector P*''^ of 
the scattered beam 

w + MP'"'' 

P- ^ (aU^Trip-a)^ c + v P^nC ^^a) 



where CTtot is the total scattering cross section. 

As can be seen, all physical quantities are expressed 
in terms of the coefficients Uki^p) which define the scat- 
tering amplitude i{(p) in Eq. To calculate these un- 
known coefficients Uki^p) for our scattering problem we 
apply the method of partial waves, similarly as in Ref.m 
Choosing the spin quantization axis along the z axis, the 
eigenspinors (the eigenvectors of the Pauli matrix az) are 
7ct, where 7+ = (1,0) for a = +1, and 7_ = (0, 1) for 
a = —1 (T stands for the transposed of vectors) . Here- 
after, we write ± for the spin quantum number a ~ ±1. 
Since the Hamiltonian Ti. in Eq. ^ commutes with the 
total angular momentum operator Jz — —ihd^ + ^az, 
any partial wave, which is a solution of the Schrodinger 
equation can be labeled by the quantum number j € JJ 
and the spin quantum number a of the incident electron. 
Here J = {• • • , — |, — ^, |, |, • • • }. We chose the direc- 
tion of the propagation of the incoming plane wave along 
the X direction. Then, in polar coordinates, the incoming 
plane wave with spin quantum number a and energy E 
can be expanded in terms of partial wavesiSi: 



,(r) = e"=--7. = 
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(8) 

are the outgo- 



where hf^'^\Y) = i/ji'^yfcr)e^(j-'"/2)^7 
ing (superscript 1) and incoming (superscript 2) waves, 

and A: = |k| = \/2m*E /h is the magnitude of the wave 
f 1 2) 

vector. Here Hm {z) are the 1st and 2nd kind of Hankel 
functions of order m. 

First, we consider the individual partial waves. The 
partial waves outside the scattering region (r > a) have 



the form i^f} = hfl + si'^hfl + S^ll^.hfl^, while in- 
side the scattering region (r < a) they can be written as 



(R) _ /lO) 



A+[aX3,+ + A%X3,-, where 



Xj,r(r) 



J,+i/2(grr)e^'^/2 



(9) 
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and Jm{x) is the Bessel function, q^- = y^W^Ht^— rfcso, ie for P'"'^ and 



aim* /h and r 



The coefficients A 



±1 is the spin branch inde3*i&. 
can be calculated from 



From Eq. ifTT)) it follows that 



and S^l]^ 



the boundary conditionsii: 
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= {dr - ikso<Jip) -00^'' 

1 r—a 


(10b) 
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valid for all j G J and a — ±1, and a^p — ^snupux + 
cos (/3 Gy . For a given j these equations involving Bessel 
and Hankel functions result in eight linear inhomoge- 



neous equations for the eight unknown coefficients A. 



and S. 



(i) 



The explicit forms of the equations are in- 



dependent of the angle (p. These coefficients can easily 
be calculated numerically. Note that the following ex- 



act relations hold S, 



(J) 



S 



(-J) 

— (J, — O" 



and S'_ 



0) 



s. 



s 



for all j G J and cr = ±1. To proceed 

further, we assume that the coefficients 5^: _(_ are known 
from numerical calculations. 

Outside the scattering region the complete wave func- 
tion describing the scattering of the plane wave can be 



decomposed as ip^^ 

, (sc) 



J6J 



(N) _ 



(sc) 



where 



t/Jo- is the scattered wave. It is easy to show that 



■>p^^ and -tp^' 



(R) _ 



(R) 



satisfy the boundary condi- 



tions H10|l . Using the Hankel's asymptotic expansions'^ 



we have /ij^^(r) 



_2_ j(fer-(j-<T/2)f ) t{j-a/2)<f 



irrkr 



To- 



valid for r ^ a, and then the asymptotic form of the 
scattered waves il>a'^^ can be calculated. Finally, Eq. Q 
yields ij^l^i^r^) = ^ 



, and using Eq. Q we obtain 



Mf) =^Bj cos (j + Cj cos (j + (^,(lla) 



ui{^) = 2sin((^/2) ^ i?,cos(jV), (lib) 
M^) = -'2cos{ip/2)^ DjCos{j(p), (11c) 



where we introduced the notations Bj ~ {S^^_^_~l), Cj ~ 

(5i^',l - 1), - and J+ = {1, |, . . . }. 

In numerical calculations the two dimensionless pa- 
rameters characterizing the scattering process are ka and 
ksoO" Figure El shows the asymmetric (skew-scattering) 
feature of the differential scattering cross section calcu- 
lated from jsj for different spin polarizations P'"'^ of 
the incident electron beam. The scattering cross sec- 
tions at a given scattering angle (p are different for spin 
up and spin down polarization of an incident beam. 
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FIG. 2: Differential scattering cross sections (in units of a) as 
functions of the scattering angle if (in units of degrees) for spin 
polarizations P'^^^ = 0, (±1,0,0)^ (a), P''"' = 0,(0, ±1,0)^ 
(b) and P'"'^ = 0, (0, 0, ±1)^ (c) with dotted, solid and daslied 
lines, respectively in each figures. The parameters are ka — 1 
and fcsocs = 1. 

u{ip = 0) ~ (0, 1,0)"^, therefore the optical theorem {7)) 
implies that the polarization dependence of the total scat- 
tering cross section (the areas under the curves in Fig.|21) 
takes the form crtot(P''"') = ci -^-C2P^"^ where PI'"' is the 
y component of the polarization vector of the incident 
wave, and ci and C2 depend only on ka and fcgoa. Thus, 
in Fig. and c, and for P'"'^ — the total scattering 
cross sections are the same. 

The magnitude of the polarization |P'''^| of the scat- 
tered waves obtained from Eq. © and the differential 
scattering cross section are plotted in Fig. |31 for unpo- 
larized (P'""^ = 0) incident electron beam using experi- 
mentally relevant parameters^S. As can be seen from the 
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FIG. 3: The magnitude of polarization IP^'^] (main panei) 
and the differential scattering cross section da /dip (left inset) 
as functions of the scattering angle (p (in units of degrees) for 
unpolarized (P'"'^ = 0) incident waves. The large scafe p (in 
units of degrees) dependence of IP'^'^I is shown in the right 
inset. The parameters are ka — 20 and ksoO, — 4. 

figure, the differential cross section is rather high and 
the scattered beam is almost fully polarized for a narrow 
window of angles. This result suggests an effective tool 
for spin-filtering without using magnetic field. 

We also studied the high energy limit of the scattering. 
This is the case, when /ca 3> 1, ie the Fermi wave length of 
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the electron is much smaller than a. Then one can apply 
the first Born approximation and it can be shown that 
the scattering amplitude f ^ ay. Thus, from ^ we have 
Mo = and u ~ (0, 1, 0)^, while from © w = 0. There- 
fore, for unpolarized incident electron beam (P'"'^ = 0) 
the polarization of the scattered waves is negligible in the 
high energy limit. However, a finite spin polarization can 
arise even in the first Born approximation if an additional 
electrostatic potential is present beside the spin-orbit in- 
teraction, eg, in an imaging experiments^. In this case. 
Mo and M2 are finite and it yields a finite value of P*"^ for 
unpolarized incident beams. 

In the opposite limit, ie for low energy limit (fca <C 
1) it is also possible to derive an analytical result for 
the scattering amplitude and the polarization P*"^ of the 
scattered waves. Keeping only the first order terms in 
k/kso of s2^^. resulting from the boundary equations l(TT)| 
and that of Uk{^p) in Eq. Hll() . it yields 

p-((p) ss2 A(_sin(^,l + cos¥>,0)^, (12) 



valid for fca ^ 1 and k <C kso and for unpolarized inci- 
dent electron beams. We found an excellent agreement 
between this result and that obtained from with nu- 
merically exact calculations. Similarly, it can be shown 
that the differential scattering cross section is approxi- 
mately isotropic (independent of the scattering angle ip) 
in the low energy limit. 

In conclusions, we have shown that for a circular shape 
of region with non-zero Rashba coupling strength, the 
scattering properties are strongly anisotropic, and the 
scattered wave is well polarized in a narrow window of 
the scattering directions. Such a nonuniform SOI can 
be utilized to produce spin-polarized electrons in an all- 
electric realization of spintronic devices. 
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